#################################################################
##DO SELF-REPORTING REGIMES MATTER? EVIDENCE FROM THE CAT########
## Cosette D. Creamer and Beth A. Simmons #######################
## February 2019 ################################################
### Replication Code for Appendix C & Appendix D ################
#################################################################
#################################################################
## Clear environment
rm(list=ls())
sink("CAT_AppendixCD")
### Libraries
library(ggplot2)

##############################################
### Figure C1: Domestic CSO Shadow Reporting #
### by CIRI Torture Score ####################
###############################################
c1 <- read.csv("figureC1.csv", header=TRUE)
postscript("FigureC1.eps", width=5.5, height=3.5, horizontal = FALSE, 
           onefile = FALSE)
boxplot(c1$shadowdom ~ c1$tort, 
        xlab = "CIRI Torture Score", ylab = "Shadow Reports (Domestic CSOs)",
)
dev.off()

##############################################
### Figure D1: Domestic Media Coverage  ######
### in Mexico and Venezuela ##################
##############################################
d1 <- read.csv("figureD1.csv", header=TRUE)
figD1 <- ggplot(d1, aes(x=Years, y=Value, fill=Reference))
figD1 <- figD1 + geom_bar(stat="identity") + theme_bw() + scale_fill_grey()
figD1 <- figD1 + facet_wrap(~Country, ncol=1)
figD1 <- figD1 + xlim(-5, 5) + labs(x = "Years Since Review", y= "Average Number of Articles")
figD1
ggsave("FigureD1.eps", width = 6, height = 5, units = "in")

################################################
### Figure D2: Domestic Media Coverage in ######
### Argentina, Brazil, Chile, Colombia, ########
####Nicaragua, and Peru ########################
################################################

d2 <- read.csv("figureD2.csv", header=TRUE)
figD2 <- ggplot(d2, aes(x=Years, y=Value, fill=Reference))
figD2 <- figD2 + geom_bar(stat="identity") + theme_bw() + scale_fill_grey()
figD2 <- figD2 + facet_wrap(~Country, ncol=2)
figD2 <- figD2 + xlim(-5, 5) + labs(x = "Years Since Review", y= "Average Number of Articles")
figD2
ggsave("FigureD2.eps", width = 6, height = 5, units = "in")

sink(NULL)
